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ABSTRACT 


The  cross  section  for  the  simultaneous  capture  of  two  electrons  by  fast 
alpha  particles  in  helium  is  computed  in  the  first  Born  approximation,, 
Calculations  are  performed  for  capture  into  the  ground  state  (  IS  )  and  the 
31SZS  state  of  helium,  In  both  cases  the  complete  interaction  Hamiltonian 
is  retained  and  it  is  seen  that  at  the  lower  energies  and  for  ground  state 
capture  the  matrix,  elements  corresponding  to  the  nucleus-nucleus  interac¬ 
tion  are  comparable  in  magnitude  to  those  of  the  Brinkman-Kramer s  type„ 
With  the  aid  of  the  Feynman  technique  are  obtained  integral  expressions, 
suitable  for  numerical  evaluation,  for  the  matrix  elements  describing  cap- 
ture  into  the  IS2,  31S2S,  31S2P^  ,  hsZP^.  ,  31S2P&  ,  and  11S2P2. 

states,,  By  the  subsequent  use  of  Gaussian  quadrature  integration  the  matrix 
elements  for  IS  and  1S2S  capture  are  evaluated  numerically  on  an  IBM- 
1620  digital  computer,, 
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Aside  from  the  inclusion  of  the  1S2S  state  in  the  capture  cross  section 
this  work  differs  from  previous  calculations  in  that  a  two-parameter  wave 
function  is  used  for  the  ground  state  of  helium »  Unfortunately  the  existing 
experimental  data  provides  only'  a  few  measurements  of  the  total  cross 
section  under  consideration  and  at  energies  where  the  first  Born  approxi¬ 
mation  just  begins  to  be  applicable,, 
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“I"  Introduction. 


The  first  evidence  of  charge  exchange  collisions  by  high  speed  particles 
passing  through  matter  was  obtained  by  Henderson*,  who  observed  that  a 
beam  of  alpha  particles  contains  a.  fraction  of  singly  charged  particles  which 
increases  with  decreasing  particle  velocity,  Since  then  a  wide  variety  of 
experimental  and  theoretical  research  has  been  carried  out  on  the  general 
phenomena,  of  charge  changing  collisions,  .A  good  summary  of  the  expert- 
mental  results  until  1952  is  presented  by  Massey  and  Burhop  . 

Due  to  the  inherent  complexity  of  electron  capture  problems  (  i,  e,  as 
pointed  out  by  Bohr  they  are  three-body  problems  in  the  simplest  cases  ) 
the  theoretical  considerations  have  b©en  rather  limited  in  scope.  Although 
several  attempts  were  made  to  calculate  electron  capture  cross  sections 
from  a  classical  point  of  view,  it  was  not  until  1928  that  Oppenbeimer^  and 
especially  in  1930  when  Brinkman  and  Kramers”  (  hereafter  referred  to  as 
BK  )  reconsidered  the  problem  quantum  mechanically.  Following  Oppenheimer , 
BK  only  considered  the  (  incident  particle  )  -  (  atomic  electron  )  interaction 
as  the  perturbation  causing  the  charge  exchange.  The  cross  sections  obtain¬ 
ed  by  this  method  were  approximately  five  to  ten  times  greater  than  the  ex¬ 
perimental  values  of  that  time,  Later,  in  1.953,,  Jackson  and  Schiff^  (  here- 
after  referred  to  as  JS  )  and  Bates  and  Dalgarn©  (  hereafter  referred  to  as 
BD  )  independently 5  and  by  different  techniques,  obtained  results  in  excellent 
agreement:  with  experiment  for  the  problem  of  single  electron  capture  by 
protons  in  atomic  hydrogen.  In  contrast  to  the  method  of  BK,  BD  and  JS 
retained  the  total  interaction  Hamiltonian  as  the  perturbation  responsible  for 
the  charge  exchange,  G,C,  Wick,  quoted  in  a  footnote  by  JS.,  argued  that 
the  (  incident  particle  )  =  (  atomic  nucleus  )  interaction  should  yield  a.  neg- 
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1  i gib  1  e  contribution  to  the  cross  section  since  in  the  limit  M/m  1  the 
motion  of  the  heavy  ions  can  be  treated  as  classical  straight  line  trajectories. 
In  other  words  this  interaction  is  a  known  function  of  time  and  can  be  factored 
out  of  the  wave  equation  by  an  appropriate  canonical  transformation.  However 

O 

Jackson  advanced  a.  counter  argument  in  which  he  demonstrated  that  at 
energies  where  1,  the  transverse  momentum  impulse  given  to  the  atomic 

electron  is  comparable  to  the  transverse  momentum  impulse  due  to  the  Cou¬ 
lomb  interaction  between  the  two  colliding  ions.  Hence  one  might  expect 
that  near  these  energies  the  Fourier  components  of  the  separate  interactions 
would  interfere  destructively  thereby  lowering  the  cross  section,  Jackson 
also  demonstrated  that  the  reason  for  the  JS  and  BD  calculations  agreeing 
so  well  with  experiment  was  due  to  the  fact  that  when  the  nuclear  charges 
are  near  unity  and  the  complete  interaction  Hamiltonian  is  retained  then  the 
second  Born  approximation  term  nearly  vanishes. 

One  serious  objection  to  the  claims  of  the  BD  and  JS  results  is  that  they 
were  compared  to  the  experimental  cross  sections  for  protons  in  molecular 
hydrogen.  Furthermore  their  cross  sections  do  not  seem  to  have  the  correct 
asymtotic  behavior  in  the  limit  of  high  incident  particle  energies.  As  a 
result  attempts  have  been  recently  made  by  McCarrol^,  McDowell^,  and 
others  to  include  the  effects  of  distortion  produced  on  the  scattered  waves 
due  to  the  (  atomic  nucleus  )  -  (  electron  )  potential.  These  results,  although 
differing  from  each  other  slightly  in  magnitude,  all  seem  to  approach  the 
BK  cross  section  at  high  energies.  However  it  will  not  be  until  experiments 
are  performed  on  atomic  hydrogen  (  if  possible  )  that  light  will  be  shed  on 
the  merits  of  any  one  approximation. 

Contrary  to  the  difficulty  implied  above  of  performing  charge  exchange 
experiments  in  atomic  hydrogen,  there  is  no  such  problem  when  considering 
atomic  helium  as  the  target  material.  Thus  in  this  thesis  the  Born  approx- 


imation.  as  applied  by  TS  will  be  used  to  compute  cross  sections  lor  the 
simultaneous  captur  e  of  two  electrons  by  fast  alpha  particles  in  helium. 

The  implication  is  not  that  the  experimental  work  in  this  regard  is  complete 
but  that  it  might  be,  perhaps,  easier  to  perform  than  for  the  case  of  protons 
in  atomic  hydrogen.  In  fact  only  a.  few  experimental  points  have  been  obtain- 
ed,  by  Allison  ,  and  these  lie  in  the  region  where  the  first  Born  approxima¬ 
tion  just  begins  to  be  valid. 

The  first  attempt  to  compute  cross  sections  for  the  simultaneous  capture 

1  2 

of  two  electrons  by  alpha  particles  in  helium  was  by  Betts  who,  after  em¬ 
ploying  the  formal  method  of  JS,  evaluated  the  rather  cumbersome  matrix 

|  'l 

elements  by  various  approximations.  Later  Gerasimenko  and  Rosentsveig 
performed  a  numerical  computation  for  capture  into  the  ground  state  of 
helium  by  approximating,  as  did  Betts,  the  ground  state  with  a  simple  one- 
parameter  hydrogenie  wave  function.  The  cross  sections  obtained  were  an 
order  of  magnitude  smaller  than  those  of  Betts,  Here  we  investigate,  again 
numerically,  the  effect  on  the  cross  section  by  employing  a  slightly  more 

complicated  wave  function  for  the  ground  state  of  helium.  The  contribution 

*2 

to  the  cross  section  from  the  J1S23  state  is  also  computed  and  the  matrix 
elements  for  capture  into  the  triplet  and  singlet  1S2P  states  are  indicated. 
Finally  a  few  possible  refinements  on  the  method  are  discussed. 
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-II-  Rearrangement  Collisions 

The  formal  theory  of  rearrangement  collisions  ha.s  been  treated  in  an 

14 

elegant  and  precise  operator  formalism  by  Lippmann  Here  we  will  review 

8 

the  summary  which  Jackson  gave  of  .Lippmann” @  paper  and  then,  apply  the 
re  suits  s  in  the  accuracy  of  the  first  Born  approximation,  to  rearrangement 
collisions  of  alpha,  particles  with  helium  atoms. 

It  is  customary,  in  a  formal  treatment  of  rearrangement  collisions,  to 
spea.k.  of  an  "initial"  and  a  "rearranged”  system 0  The  Hamiltonian  for  the 
total  system  may  be  written  in  either  of  two  ways: 

H  -  Ho  +  v/  =  H J  V1  (2-1) 

where  the  unprimed  quantities  refer  to  the  initial  system,  the  primed 
quantities  refer  to  the  rearranged  system,  HQ  and  Hc"  are  the  respective 
"unperturbed”  Hamiltonians,  while  V  and  V"  are  the  respective  "inter- 
action’  parts„  The  unperturbed  Hamiltonians  have  eigenvectors  <§co  and 
such  that 

Ho  $(X=  j  9b 

The  operator  Green’s  functions  Gc  and  G0”  are  defined  by 

&o(  x'l  -  vc~h7  ,  Go'<  ^  H.' 

where  X=£.±Le  for  outgoing  (  +  )  or  incoming  (  waves.  The  solutions 
describing  incident  waves  plus  scattered  waves  satisfy  the  integral  equation 

=  ?«,  +  (2-4) 

and  the  corresponding  primed  equation.,  If  one  defines  an  operator  U  (  AN) 
such  that 


(2  =  2) 


(2=3) 


l  1  -  IKV)  =  1 


(2-5) 
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then  equation  (2-4)  (or  its  primed  equivalent)may  be  written  more  compactly 
as 

(±) 

fee  =  U(Eatt£)g>a  .  (2-6) 

Lippmann  shows,  rigorously,  that  the  transition  probability  for  a  rearrange- 
ment  collision  is 

c*V  cw)  (2-7) 

and  furthermore,  that  if  both  and  H0’  are  Hermitian 

=  <tr»vnt>  ;  £b'sE^  (2-8> 

and  also 

<  fb'  IV\  So.)  ~  <ft'lvl&u>  •  El-Zc  (2=9) 

^■r.  C+“ J> 

The  first  Born  approximation  consists  of  replacing  '£0v  by  in 

eqn„  (2  =  7)„  This  corresponds,  formally,  to  a  neglect  of  all  terms  but  the 
first  in  the  expansion 


U  (  ^  "  1  h  C~0l  XW  “  G0(  vw  Ge(>.'w  +  *  '  -  (2=io) 

or  its  primed  equivalent0  Thus  the  differential  cross  section  for  a  rearrange¬ 
ment  collision  in  which  the  initial  configuration  is  characterized  by  and 

the  rearranged  system  by  ?w  is  given  by  the  familiar  formula 

C  ©»  $) s  )  vf-  '^bcJ  (2-ii) 


where  is  given  by  either  expression  in  equation  (2  =  9) „ 

To  specialize  to  the  specific  problem  of  two  =  electron  capture  by  alpha. 

particles  in  helium  we  utilize  the  coordinate  system,  illustrated  in  Fig0  1, 

1  7 

used  by  Betts1  in  the  same  connection  The  nuclei  are  designated  by  A  and 
B,  1  and  2  are  electrons  and  and  C  are  the  centers  of  mass  of  the  initial 


and  final  atoms. 


■ 

' 


’  'I 


•;  i  'i 
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Fig0  1  Coordinate  System 


If  we  consider  ion  B  initially  incident  with  relative  wave  vector  upon 

the  two-electron  atomic  arrangement  about  A  and  finally  a  two-electron 
arrangement  about  B  departing  from  ion  A  with  relative  wave  vector  , 

we  may  write 


where  and  X|lfBi)rsxN)  are  the  initial  and  final  helium  atom  wave 

functions  respectively. 

Making  the  usual  assumptions  that  11?^!”!^.!  and  •&  5  the  scattering 
angle,  is  small  the  momentum  transfer  may  be  written  as 


*  Ru  ~  *  R/ 


rBl.  A  +  r^x>  (a  tfi)  -  B 


where 


OU  A  »  (.VT I1'  s1 }  =  4  i  r0  )  (“Ar) A  •  B  =  O  (2-14) 

and  is  the  velocity  of  an  electron  in  the  first  Bohr  orbit  of  the  hydrogen 


!  ■■  .■  v'J.isi-i.  ;  :rr  cl  ■  tabia  ;>  -•«!  >] 


"-■  ■  r  " ■  ' ">■  •  i'ii? p  it*-?  x.&  ...rcnoii  nc<"i. r» -  3>/;?  -a-rfi 

"  •  >  it ■■■:  "  ■;  :  i  8  7(.'C  C  -•  i.r.7  ?.v  .  7  ' 


X  . 


3'.  •  ,  '  '  A  ” 

SI 

:  1  ■-  ••  >s  ■  .7  '  •  ■■ ;  ;  ,  ,/  -c  ■; 


"  uf-J 


3  ff 


1  ^  iao£  ir  .  ■  ij  a  >••  ,;a  **«  ‘  ,  yti:  -.(A..  .  v  biw 
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atom,  Thus,,  dropping  the  subscripts  b  and  a„  we  have 


X  £xpt[rB,*K  +  rM-(^+B)-  rSl.g]  d?A2  3>s,'3*' 


Since  V  (.  R  >  ~  ~JL§t  4-  » 

rA*  R 

write  (2-15)  more  conveniently  as  the  sum  of  two  terms;  thus 


Bi  4  (2= 15) 

,  we  may 


^  “  H^bk)  V  H^moc^  (2=16) 

where  is  the  Brinkman=Kramer s  type  matrix  element  given  by 

Hu*')  =  -Ae1/ fy*(rBl)FeO  X,4(“AI)7„\ 

X  exp  +^(f,+‘6')- vC_L_  &A1  3r»,  3vSi  (2- 17) 

rh\ 

and 

^  (.NOt')  “  (.  'Oj,*  V'^z')  l  rA*  > 

X  4  (At  s')  -  %,;  b1  V-  t 

a 


. 


- 
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III  Some  Simple  Hydrogenic  Wave  Functions  for 
the  Ground  State  and  the  First  Few  Excited 
States  of  Helium 


In  this  section  we  shall  review  some  familiar  results  from  the  variational. 

method  of  constructing  approximate  eigenfunctions  for  the  helium  atom . 

As  well  we  shall  present  the  results  of  some  calculations  of  'effective  charge' 

3  3 

parameters  used  in  simple  variational  approximations  to  the'  1S2S,  '  1S2P, 
and  ^  1S2P  helium  wavefunctions0 

The  well-known  variational  method  consists  of  first  constructing  a  trial 

function  of  an  arbitrary  number  of  adjustable  parameters  ,  say  , 

which  satisfies  the  same  boundary  or  asymtotic  conditions  as  the  true  eigen- 

function,.  If  T|f  is  normalized  to  unity  and  constructed  so  that  it  is  orthogonal 

1 5 

to  the  first  j-1  states,  then  it  can  easily  be  shown  that 

"jVV  <^!H  I  f)  ^  Ej  (3=1) 

where  Ej  is  the  true  energy  of  the  j  th  level  and  H  is  the  total  Hamiltonian 
of  the  atom  in  question.  The  Hamiltonian  for  the  helium  atom  may  be  written 


as 


H  -  -TTCfo 


?!)  -zel(i  + 


r,  -  p.  I 


(3-2) 


where  r,  ,  r^are  the  coordinate  vectors  of  the  two  electrons  and  and 

.2, 

^  are  the  corresponding  Lapiacian  operators.  Since  the  inequality  in 
eqn,  (3  =  1)  is  always  valid,  provided  satisfies  the  stipulated  boundary, 
normality,  and  orthogonality  conditions,  we  say  that  the  'best®  choice  of 
the  parameters,  ^  ,  is  obtained  by  the  n  conditions  which  minimize  the 
energy,  E(<x,j*a  }  *• 


'd&C 


0 


<■'  J  a  • 


(3-3) 


3  °<, 
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The  above  Hamiltonian  is  invariant  under  the  group  consisting  of  the  two  dis¬ 
crete  operations  of  permutation  of  the  two  electrons  and  inversion  of  the  coordi¬ 
nate  system  and  the  continuous  set  of  operations  of  rotation  of  the  coordinate 
system,  As  shown  for  example  in  Wigner^®  the  eigenfunctions  of  the  Hamiltonian 
act  as  basis  functions  for  the  various  irreducible  representations  of  the  group  of 
H„  Hence  even  if  we  cannot,  a.s  in  this  case,  find  the  true  eigenfunctions  we  do 
know  what  symmetries  they  possess.  Further  we  are  guaranteed  by  group  theo¬ 
retic  considerations  that  functions  of  different  symmetry  in  the  sense  that  they 
belong  to  different  representations  of  H  (or  different  rows  of  the  same  represen¬ 
tation  in  the  case  of  degenerate  eigenvalues)  are  automatically  orthogonal.  Thus 
the  variational  method  can  be  used  to  find  an  upper  bound  to  the  smallest  true 
eigenvalue  corresponding  to  each  different  symmetry. 

Because  of  the  limitations  of  time  and  because  as  in  the  protons  in  hydrogen 
case  only  the  lowest  states  are  expected  to  be  important,  we  consider  only  the 
subgroup  of  the  discrete  operations  of  inversion  and  permutation.  Under  inver¬ 
sion  the  wave  functions  may  be  invariant  (even  parity)  or  change  sign  (odd  parity) 
and  likewise  under  interchange  of  coordinates  of  the  electrons  the  wave  functions 
may  be  invariant  (symmetric)  or  change  sign  (antisymmetric).  In  atomic  physics 
the  spectroscopic  notation  for  the  lowest  states  of  each  of  these  four  possible 
symmetries  are 

IS  ‘  even  symmetric, 

3 

1S2S  even  antisymmetric, 

^TS2P  odd  symmetric, 

3 

1S2P  odd  antisymmetric  „ 

The  only  other  state  of  comparably  low  energy  is  the  ^  1S2S  state.  However  as  it 
is  of  the  same  symmetry  as  the  unknown  IS  state  the  variational,  method  cannot 
guarantee  an  upper  bound  to  the  ^  1S2S  energy  level  and  hence  it  is  not  possible  to 
say  how  good  the  variational  energy  levels  and  wave  functions  for  this  state  are 
and  so  we  ignore  it  in  the  remainder  of  this  thesis. 
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A„  The  ground  state  (IS1-); 

The  simplest  hydrogenic  trial  function  for  the  ground  state  of  helium,  as 
shown  in  almost  any  text  on  quantum  mechanics,  is 

s/x.  -2'/q0(i",  +  >\) 

ijr  -  (  ^/rrQfl  )  6  (3-4) 

which  when  substituted  into  eqn„  (34)  yields 

EC  2')  =  Inc  „  (3-5) 

The  subsequent  minimization  of  E(g)  leads  to  the  well  known,  result  that 

E,-  -  5.WS3  i?  »e  Vic.  tor  2'=  a7/ \<o  . 

A  slightly  more  elaborate  trial  function  for  the  ground  state,  in  the 
notation  of  Pauling  and  Wilson  ,  is 

(3-6) 

1 7 

where  N0(ZjH'’)  is  a  normalization  factor,,  Eckart  and  Hylleraas  used  this 
function  in  a  variational  treatment  and  obtained 

E,=  -5.1So6  }  ■2"=-  \.iq 

as  compared  to  the  observed  energy  of  -5*8074  R_Hc  for  the  ground  state  of 
helium,, 


B„  The  ^1S2S  state; 


Here  we  choose  as  a  trial  function 


^T^lSzs  ”  "^$2" Cl)  1^2*C^^3  (3  =  7) 

where  N  fe'g’*),  the  normalizing  factor,  is  given  by 

i  v -  V* 

Hu'>2c')  -  7^**  (  1  -  K) 


and 


(Z,Z,|f^(Z/°  2") 

A~  zfH  (2»  +  &/z  » 

Substitution  of  (3  =  7)  into  the  energy  integral  (3-4)  yields 


(3  =  8) 


-11 


C*-Ax) 


(2’%  E"7^y  ~U*  '  +  &'•)  +  ZAS(4“2'-  2")  +  P  ~  £ 


(3-9) 


where 


B  =  n/a  (Z'Z'THz^  I?") 


(3-H 


Pc  £Q./lS*(l);LS*U}  lFt-\\  lSaiCl)aSa*U5d^,  ^ 


A. 

u 


Cz'2")  /  4  .  ^  4 \ 

2»+2*^5  +  2DZIZ1'  +IH1  Z"  +10Z'2‘*  y  (3-11] 


and 


Q  -  xao /  iv  (i)  2^,.^  jTTTJ  s-s^ul  &?, 


_  (z'A  [oo.  _  a-PQ  *"  ,  13a.  •gl'a'  1 

IfoU'  +  ’gj'f  L  U 1  +  t!  xj  '  (z»  +  l."  y  J  ,  (3-12) 

The  reader  is  referred  to  Appendix  1  for  the  method  of  evaluating  the 
integrals  P  and  Q, 

Obviously  an  application  of  the  minimization  conditions  (3-3)  to  eqn„ 

(3  =  9),  would  be  a  formidable  task,  Instead  we  argue  that  it  is  physically 
reasonable  to  expect  that  eqn,  (3  =  9)  would  be  minimized  for  ?•'  somewhere 
in  a  region  about  =  and  in  a  region  about  Z'1  - 1  „  With  this  in 

mind  a  digital  computer*  was  programmed  to  evaluate  eqn„  (3  =  9)  for  a  wide 
range  of  points  with  2 1  and  being  incremented  (or  decremented)  in 

steps  of  0,  01„  Eventually  the  deepest  minimum  was  found  yielding  the 
values 

E«  -  4.3333  Rh«,  he  5  *'*  l.<W  >  S 

as  compared  to  the  experimental  energy  of  ^  for  the  ^1S2S  state 

of  helium, 

- - — — - — - n- - - - — - - - - — — 

*  The  University  of  Alberta's  IBM~1620j  and  LGP-30  were  used  for  the 
variational  calculations  of  this  chapter. 
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C The  ^1S2P  and  ^IS2P  states  ° 


Consistent  with  the  form  of  the  previous  trial  functions  we  write  here 


Substitution  of  eqn„ 


(3-13)  into  (3-1)  yields 


(3-13) 


=  1(2'"*  *7+)  -Mi'+e')  +  r  +  s]RHe^ 


C  (3- 


where 


R  ~  Pe»'U)  -  rj  ^  p2„(.2-\<Jr 


and 


_  Z^CSZ'^  +  IR-Z'V  +IS'Z,E"’'+  st'?) 
&i-A  l  ZH  +  10s" 


(3-15) 


ISgiU.)  i  p^tf)  ^7 


_  _7_  7'a&"r 

t*  (?'  +-  4."V  •  (3-16) 

Again  the  reader  is  referred  to  Appendix  1  for  the  method  of  evaluating 
the  ’coulomb”  and  ’exchange"  integrals  R  and  S  respectively„ 

Finally,  numerical  minimization  of  eqnQ  (3=14)  yields  the  results 

E.Msip  z‘a\.^8  , ,  z"»  l.os 

EMSIP  *  -  4-5.443  RHe>c  ^  ^  -£.oo  )  z"  =  i.oo 

compared  to  the  experimental  values  of  ~  4‘2.(eb3  Rv^gloc  and 
""  4*5-47  b  KC  respectively,, 
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-IV  ■  Cross  Sections  for  the  Simultaneous  Capture 
of  Two  Electrons  Into  the  Ground  State  by 
Fast  Alpha.  Particles  in  Helium 

1  2 

The  simple  wave  function,  eqn0  (3-4),  was  first  used  by  Bells  and  later 

1 3 

by  Gerasimenko  and  Rosentsveig  (hereafter  referred  to  a,s  GR)  in  a  computa¬ 
tion  of  cross  sections  for  two- electron  capture  into  the  ground  state  of  helium „ 

The  cross  sections  obtained  by  Betts,  based  on  an  approximate  evaluation  of 
the  BK  and  nuclear  matrix  elements,  are  an  order  of  magnitude  higher  than 
those  obtained  by  the  semi-numerical  procedure  of  GR  „ 

To  serve  as  a  ’starter'1  for  the  more  elaborate  calculations  which  follow  we 
have  also  made  a  numerical  computation  of  ground  state  capture  cross  sections 
using  the  same  wave  function  as  Betts  and  GR  „  The  pr e-numerical  manipulations 
lead  to  matrix  elements  expressible  as  integrals  over  a  unit  cube  whereas  those 
of  GR  are  over  a  unit  square.  However  the  relative  simplicity  of  our  integrands 
compared  to  those  of  GR  is  sufficient  justification  for  the  extra  dimension, 
especially  in  view  of  the  more  complicated  wave  functions  which  will  be  used 
later  on. 

The  procedure  is  as  follows;  inserting  expression  (3=4)  into  eqn„(2=17) 
we  have 

HCbi0=  "^7  exp[-°<(rM +rA1 +rsi -i-r62.) +r„-(A  +  %’)-f82-s)i 

X  ~rZ  3rBl  (4-1) 

where  c<-  ^I\(qQ0  »  Introducing  the  Fourier  transforms  (viz.  Appendix  2) 


14 


e 


-  oL  r 


A1 


°L  (  e  til '  do 


(o^  +-Q,Z 


etc . 


into  eqn0  (4 -l)  gives  us 


V 


H 


L<6<) 


2_  cX 


IT' 


€xp  6(  rAz.  (^|+A +B+^)  *  (|a  +  A )  +  rBz  *  (^>  6  -  ^  )  J 


— 1 

fv^ei 

© 

o^+f  'l 

1  ( f  U* + y  7  0*1 + ?f  )z 

(4-2) 

Making  use  of  the  Fourier  transform  of  the  cT  ~  function  8 


dU  r  -  r0)  = _ 1 

eqn„  (4  =  2)  readily  reduces  to 

^(5^  [ - 


c^y 


®  (  r  to  )  ^ 


T-  J  +  £■  )  [o< 1  4.  (  ^  +  A  |>  h  t|4  a)^[  +  (  o  +  -%fj 


1  8 

By  the  use  of  a  trick  introduced  by  Feynman  in  another  connection  we  can 
reduce  eqn„  (4  =  3)  to  a  form  more  appropriate  to  numerical  computation,, 

Each  of  the  factors  in  the  denominator  of  (4=3)  is  expressed  in  terms  of  symbols 
a9  b5  c,  d  and  use  made  of  the  identity  (vize  Appendix  3  ) 

t  />  i  /'i 


(X  bx  Cz 

where  we  put 


izof  f  I  Y^Ci-xXl  -  yti-  dxdycU _ 

[a^cyz  -I-  byz(v-*)  f  Cz  (\~y)  J7 


a.  =. 

cX2- 

+-  c 

b  = 

h 

(|  h  A  +b) 

c  = 

of2- 

4 

(f  4Af  , 

ii 

—n 

o<Z 

4- 

-v if 

and 

Simple  algebraic  manipulations  lead  us  to  the  result 

_  sJ°x 


H 


^k) 


Tr  t72'°'> [J  I 

Jj(p  +/\f+-  ^ 


7 
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where 

A  =  AC,  -I-  e>Cx 

=  °C-  +  AaC,  (l-C,)  +  &c1_(\-c2) 

C,  =  2  (  l-Xy) 

and 

C-,Z  =  1-7  +■  7?  -  K'y? 

Finally  when  use  is  made  of  the  well-known  result 


Ia 

eqn.  (4-4)  reduces  to 


TTX  j2$-5)  U 

JJ-iL 


tk z+fiz)x  tt-iV.  5 


H  (6^-'  -  30,2.40 oc'r 


f  f  f  vx(  i-x)(i-yK  >- z)g4  dx  dy  <Jg 
JoJ0]o[  AzC,(i-C,)  +  B^O-CO  +  X1  ] 


]"/= 


It  is  noted  that  when  JS  applied  the  above  methods  to  the  protons  in  hydrogen 
problem  the  integrals  defining  the  nuclear  and  BK  matrix  elements  could  be 
evaluated  exactly.  In  the  present  case  however  we  have  two  equivalent  integrals, 
(4  =  3)  and  (4  =  6),  neither  of  which  can  be  analytically  integrated.  Inspection 
of  (4  =  3)  shows  that  the  integrand  is  very  sharply  peaked  near  q“08  q~  =  (  A  +  B  ), 
q  =  =  A,  and  q  =  =B.  This  was  the  basis  for  the  approximate  evaluation  of  the 
integral  by  Betts.  However  it  is  very  difficult  to  find  a  numerical  procedure 
which  will  take  advantage  of  the  fact  that  the  chief  contribution  to  the  integral 
comes  from  the  above  points,  especially  when  each  evaluation  of  the  integral 
will  be  for  a  different  value  of  A  or  B  or  both.  Hence  the  superiority  of  (4  =  6) 
to  (4=3)  in  a  numerical  computation  is  clear  since  in  (4  =  6)  the  peaking  is  not 
so  evident  and  the  range  of  integration  is  drastically  reduced. 

The  nuclear  matrix  element  is  attacked  in  a  similar  manner  to  that  outlined 
above.  Substituting  (3=4)  into  eqn.  (2=18)  and  again  introducing  the  appropriate 
Fourier  transforms  for  1/R,  etc.  (see  Appendix.  2),  leads  us  to 
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d  (noc) 


cfo  c^p 

<-*-  FvZ  in  .  -K  ^  o  \2- 1 2_  1 


7 T3  7  px  C  *z  [°<z  H<|  -  p  -  5 +  (|-A)^^  +  (p  p-A--B)^]4 


Eqns.  (4*3)  and  (4*7)  are  the  same  as  Betts5  (4=3)  and  (4-4).  The  reduction,  of 
eqn0  (4  =  7)  to  a  form  suitable  for  numerical  evaluation  is  outlined  a, a  follows. 

Considering  first  the  integration  over  p~ space  only  we  have 

I  - 


4e. 


p1  |_D<i+cp-cy]z'[^  +  cp-E3i]J 

where  C  =  q  =  B  and  *D  -'q  =  A  =  B„  Introducing  the  Feynman  auxiliary  integral 

\  _  f 8  fc  x  ( t  -  x  ^  dx 


o  by- 

one  obtains  by  choosing 


and 


that 


[ax  +-  bO- 

a  =  c<2  +  p*  +  C  -  a-p  ■  C 

t>  ■  cT  +•  pz  +-  o  -  2.  p  •  X) 

Te 


l"  I 


where 


and 


Sinci 


pi  P  +• 

A  =  (7  x  -t-  -  O 

f  -  xC  +■  t\-*) "5 


4 


P 


IT 


we  have  that 


p’-lpVA-lpl'y-  CA-f2)'^  A 


I 


*  (-'"X)  BA2  t  (A-11)’/2-  A  j 


=l  IT  I  X 


t>y  _iz±_  + 

(A-t^A2-  (A-1Z)‘/2^ 


+ 


d 


x  * 
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Fortunately  the  term  t  &  ~  "3  )  does  not  depend  on  the  variable  q  and  hence, 
since  the  A^S  are  raised  only  to  interger  powers,,  the  Feynman  technique 
can  be  used  again  to  treat  the  remaining  q= space  part  of  eqn0(4-7)0  We  even¬ 
tually  get 


U  -  4on  ic  f  Vf'  \ 

H^-480^JJJo  i 


(4-8) 


where 

F  =  oCz  +  *(.'-*)  A,1" 

/12  =  A*C,00  -I-  B2  C-i.C'-ti)  +- ; 

C-\-  \  -  y  +  -  my*  . 

and 


C-j.  =•  y  a 

The  quantities  c ^  and  c-,  appearing  in  eqns„  (4  =  6)  and  (4=8)  are  defined 
differently  in  each  case.  Due  to  the  enormous  number  of  definable  terms  which 
occur  in  this  type  of  work  the  policy  adopted  throughout  this  thesis  has  been  to 
assign  the  same  symbols  to  terms  which  appear  in  the  same  algebraic  manner 
throughout  the  manipulations „  Thus  following  each  final  statement  of  a  matrix 
element  there  will  usually  appear  a  list  of  symbol  definitions,. 

For  an  indication  as  to  which  numerical  procedure  would  do  the  most  justice 
to  (4  =  6)  and  hence,  it  is  presumed,  to  (4  =  8)  use  was  made  of  the  fact  that  (4  =  3) 
can  be  integrated  exactly  for  the  case  of  zero  scattering  angle  (  g=0=B)o  In 
fact 


14  IP,  \  Col  =  -  2?  TT  oC4  (35fa  olf_+  3-4  «<  S 

It  was  found  that  the  use  of  eight  points  along  each  axis  in  a  Legendre-Gauss 

19  x 

quadrature  technique  7  gave  at  least  four  figure  accuracy  in  the  range  S  =■  c< 

to  8°^  .  With  a  non=zero  g  and  a  fixed  value  of  S*  the  results  of 


z 


the  numerical  integrations  converged  to  at  least  four  figure  agreement  when 


'  .  ■■■  ■  - 
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the  number  of  abscissae  per  axis  was  greater  than  seven„  Thus,,  although  the 
nature  of  the  error  arising  from  the  application  of  the  above  method  to  the 
nuclear  integral  of  eqn0(4«8)  and  to  the  more  complicated  matrix  elements 
which  follow  is  not  quite  so  clear,  this  is  the  numerical  technique  used  from 
here  on  throughout  the  thesis  i.  e.  an  8= point  Legendre=Gauss  quadrature  for¬ 
mula.  All  matrix  elements  were  evaluated  on  the  University  of  Alberta’s 
IBM  =  1620  at  an  average  time  of  3/4  of  an  hour  for  each  evaluation  of  a  BK 
element  and  1  and  1/4  hours  for  each  evaluation  of  a  nuclear  matrix  element. 
Furthermore  each  matrix  element  was  evaluated  at  several  different  values  of 

ft 

S  and^at.  each  of  these  values  of  S  ^  over  a  suitably  wide  range  of  values 

Of  go 

The  total  capture  cross  section  into  a  given  state  is  obtained  from  eqn.  (2  =  11) 

by  making  the  approximations 

,  d -Q,  *  TT  C'M-)  d(^)  J  TSi-Vf 

and  thus 


cr 


T  Cu7- 


ir 


^Jo  i  8(8^  +  1-tcNDC.)!  • 


(4-9) 


The  methods  employed  when  using  the  improved  wave  function  (3  =  6)  to 
evaluate  the  ground  state  capture  cross  section  are  similar  to  those  abovec 

W  e  obtain 


2-  irrr 


K 


+(p6>1«f  +(^d(4=1 


where  ^ 

=  A2C,(\-C,)  + 

fi  =-  +■  a2  c.2.  C\  ~c 2.1 
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tab  - 


1 


bcx. 

c» 

c, 

a 


GA  y  2  +  b2-  C  \  ~  yz ) 

bz  y  2  4-  a.1  C.  v  -  yz) 

I  -  2  4  y^z  -  xy  2- 
i  -  yz 

<&•  1 5~  ^  b  =  l*  l  C1 


) 

; 

; 

;> 


J 


and 


M  =_ L_l  _ _ 

'4*  o  i/TzT  l  (a-yb)k  -4  C-  4  (^b)3  - 


'/i 


With  a,  b,  and  N0  defined  as  above  the  nuclear  term  yields 


H^joO)  "  4  80  NiJ  tt(gV) 

1 


f\p\  p\ 

o 


4 

K  2 

L=1 


R  '/z 

‘  t 


Lf^  +  ii)7!'- 


o  J 


3 


Jo 


chclydz  KyZ'C'-x)(\-'/X\-2) 


-t- 


14  e, 


4“ 


43  Ct 

(/5l  +  ^ 


where 


F\  =•  +  xU-x)  A2' 

Fa,  ~  03- x  +  bz(j-x)  4 

b*,  ~  bzx  +•  CAC\-x)  it  xO~x)AL 

b  4  -  CA  4-  X  C  \  ~  *")  A2- 


C >  -  u 


y  +  yz  -xy^  5  CL^y^ 


-  b^Cj.  4  a2-  Ci -co 
V  b'et  4-  a>Ci-0 
-b  tACi-O 


? 

> 

> 

*> 


and 


l-v-olC;.  4-  tKi-eJ 

Tables  1  and  2  depict  the  calculated  values  of  and  from  for~ 

mulae  (4=10)  and  (4-11)  respectively.  A  comparison  of  these  results  with  those 
using  the  simple  wave  function  obtained  above  and  those  of  GR  is  illustrated  in 
Fig.  2. 


' 
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Table _ 1_ 

The  following  is  a,  listing  of  the  values  of  times  minus  one)  as  com" 

puted  from  the  two-parameter  wave  function  for  the  ground  state  of  helium 
using  eqn„  (4--10)„ 


g 

S' 

S  *  Sod2’ 

sW 

s1-  a*2- 

0.  0 

16. 1357 

6. 7832 

3. 3494 

1. 8479 

0. 7074 

0. 3325 

0.  1 

15.4321 

6.2672 

3. 0049 

1.6141 

0. 5881 

0. 2637 

0.  2 

13. 5475 

5. 0015 

2.2211 

1. 1177 

0. 3647 

0. 1488 

0.  3 

11. 0200 

3. 5472 

1.4242 

0. 6621 

0. 1916 

0. 0715 

0.4 

8.4136 

2. 3105 

0. 8337 

0. 3584 

0.0931 

0. 0323 

0.5 

6. 1127 

1.4234 

0. 4636 

0.  1857 

0. 0440 

0. 0144 

0.  6 

4. 2813 

0. 8491 

0. 2522 

0. 0951 

0. 0209 

0. 0065 

0.7 

2. 9241 

0.4994 

0. 1369 

0. 0491 

0. 0102 

0. 0031 

0.  8 

1. 9662 

0.2935 

0. 0753 

0. 0259 

0. 0051 

0.9 

1. 3118 

0. 1740 

0. 0423 

0. 0141 

1.0 

0. 8739 

0. 1048 

0. 0244 

0. 0080 

0. 0015 

0. 0004 

1.  1 

0. 5840 

0. 0643 

0. 0046 

1.2 

0. 3929 

0. 0403 

0. 0028 

1.3 

0. 2669 

0. 0017 

1.4 

0. 1833 

0. 0011 

1.5 

0.  .1274 

0. 0113 

0. 0024 

0.0007 

0. 0001 

0. 0000 

1.6 

0. 0897 

0. 0005 

1.7 

0. 0640 

0. 0003 

1  „  8 
10  9 

2.  0  0.  0251 

2.  5  0.  0005 

3.  0  0.  0022 
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T  able  2 


The  following  is  a  listing  of  the  values  of  H^JC\  as  computed  from  the 
two  parameter  wave  function  for  the  ground  state  of  helium  using  eqn„  (4-11). 


g 

$x  1*? 

2.  t  ? 

5X-  boC- 

s*  3* 

0.  0 

11. 0890 

4.  1739 

1. 8671 

0. 9421 

0. 3074 

0.  1256 

0.  1 

10. 7458 

3. 9542 

1. 7373 

0. 8632 

0.2742 

0. 1094 

0.  2 

9. 8081 

3. 3932 

1.4236 

0. 6816 

0. 2042 

0. 0777 

0.  3 

8. 5012 

2.6993 

1. 0689 

0.4907 

0. 1386 

0. 0507 

0.4 

7. 0775 

2. 0476 

0. 7674 

0. 3399 

0. 0920 

0. 0328 

0.  5 

5. 7310 

1. 5200 

0.5453 

0. 2355 

0. 0621 

0. 0219 

0.  6 

4. 5675 

1. 1266 

0. 3921 

0. 1669 

0. 0435 

0. 0153 

0.7 

3. 6193 

0. 8450 

0. 2889 

0. 1222 

0. 0317 

0.  8 

2.  8739 

0. 6460 

0. 2189 

0. 0923 

0. 0240 

0.  9 

2. 2994 

0. 5048 

0. 1705 

0. 0719 

1.0 

1.  8601 

0.4030 

0. 1361 

0. 0574 

0. 0149 

0. 0052 

1.  1 

1. 5251 

0. 3282 

0. 0468 

1.2 

1. 2653 

0. 2719 

0. 0388 

1.3 

1. 0640 

0. 0326 

1.4 

0. 9055 

0. 0277 

1.5 

0. 7791 

0. 1677 

0. 0566 

0. 0238 

0. 0061 

0. 0021 

1.6 

0. 0206 

1.7 

0. 0179 

1.  8 

0. 0157 

1.9 

2.  0 

0. 4161 

2.  5 

0. 0526 

3.0 

0.1665 

Cross  sections  for  the  simultaneous  capture  into  the  ground  state 


of  two  electrons  by  fast  alpha  particles  in  helium:  _________ 

presents  the  improved  calculations,  ;  __  „  __  <,  __  „ 

presents  our  cross  section,  ,  arising  from  the  simple 

depicts  the  results  of  GR, 


re  = 

re  = 


wave  function  and 
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-V -■  Cross  Sections  for  the  Simultaneous  Ca.pture  of 
Two  Electrons  into  Excited  States  by  Fast,  Alpha. 
Pa.rticles  in  Helium 


The  differential  cross  section  for  a  capture  process  whereby  one  electron  is 
taken  into  an'n=l  state’  and  the  other  into  an  ’n=2  state'  is  indicated  by  the 
proportionality 

+  !  Hai52sI  51*  ilH'isiPrl  \ 

In  this  section  we  present  calculations  of  the  cross  section,  ^3\$2tS  °  Matrix 

elements  for  the  singlet  and  triplet  1S2P  capture  processes  are  derived  and  left 

in  a  form  amenable  to  numerical  evaluation , 

Substitution  of  the  wave  function  (3  =  7)  into  eqns0  (2=17)  and  (2=18)  plus  mani  = 

pulations  similar  to  those  outlined  in  the  previous  section  yields  the  following 

3 

BK  and  nuclear  matrix  elements  for  capture  into  the  1S2S  state: 


H-  lex) 


±d  l 


r-  1 

n  1 

1  ~r  i 

!!  o<z  2-  (\~y) 

I  ~ 

o  ^ 

O  - 

(/S  -  «£) 

iio 


(5  =  1) 


where 

d  \  =  -  b  =  US 

-  a.  -  z,\5  j 

D<~\-  L  S  °1  J  °<2,“  O'  7  8  y 

T  =  yH'-*)('-'/XV4)24  5 

a  --  a1gvCi-0  , 
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C;=  1-2  ;  CL-  l-^y2:  ; 

^(.p^'jT's)  s  p*y?  +  +  r^^‘7 )  +‘^C\-2r')  ; 

at=  ^a.1)  , 

=^(.a *■  >°<r  jotv1-;  *4  ,  f?4  =  a  (b’-')o<li  a1) 


and 


F  C^vl  o  c.}  =. 


n,  tss- 


6- 


/■  i  r  1  r  \ 


a  Jo  J 


F; 1/2  /*£ 


T, 


A 


_fe_  + 

FV  '  ^ 


a.8c^  +  ufcc* 


/3i 


-*tTR  (_l5_ 


where 


p ,  3 


4* 


kg>  C, 


j.  n  ,1 
o 


R"/R 


~h 


18°1  C- 


Pj  jBi* 
Z  1 


I  ^3CfUK<LJ; 


4 


WC  ( 


^fey 


i__  ^  _Ufa  . .A)  cUdy  ^ 


^6 


(5  =  2) 


Ta=  *yz(i-xYv-^)C\-f)  ^  -  *TK  ^  Tc--  (P-^  Ta  ^ 

FCtju)-  Ik  -I-  U.l\~*')  h  xlv-X^A2-  jl 

F\  ^  F  0>C ;  bO  ,  Ft-FC<**aO  ,  (5-3) 

-  F  j  t>0  5  =  F  j  cll ) 


>  J 


/S^  =  [a2c,Ci-c,')  +  BicAi-cz')  +  iPl 

C  \  -  '-7  +7?  -  X'/b-  j  C  *.  ■  y  ^ 


V, 


(,s-<0 

tar-5-) 


-26  ■ 


it  pi‘pr-s')  -  p*yf  + 

=  1  j  b*  ^  a> ;  o(v 

13  --  3  (.« ^  bs  a1 ;  o^)  ,  V-5  i<xl; b ;  t^)  , 

-J 

and  as  above 

a=  a.is  ,  b*us  j  ocv»  i.'W  ,  o<v=.  o.is  . 

Tables  3  and  4  list  the  values  of  H  and  H  calculated  from  the  above 
equations  while  Fig,  3  illustrates  a  comparison  between  the  cross  section  for 
1S^  (as  given  by  Fig,  2)  and  ^1S2S  capture  by  plotting  the  ratio  |00  ^3  \S  iS^OT^ 
against  the  energy  of  the  incident  ion,  Finally  Fig,  4  illustrates  the  cross  section 
T  t  +  CTj  (e  2<,  along  with  the  experimental  points  of  Allison, 

Since  the  1S2P  states  are  not  represented  by  spherically  symmetric  wave 
functions  we  must  specify  directions  for  the  vectors  A  and  B,  We  choose  them 
such  that  A  and  B  define  the  x  and  2  axes  respectively  of  the  final  state  wave 
function.,  Such  being  the  choice  we  have  thus  eliminated  the  contribution  to  the 
cross  section  07,2,  from  all  the  ISZP^.  states.  This  is  clear  upon  realization 
that  the  perturbations  causing  the  capture  are  rotational ly  invariant  and  that  with 
A  and  B  chosen  as  they  are  the  introduction  of  the  1S2P^  wave  functions  into 
either  of  eqns„  (2  =  17)  or  (2=18)  is  tantamount  to  averaging  a  function  of  odd  parity 
over  all  space.  Of  course  this  result  will  readily  appear  upon  a  detailed  examina= 
tion  of  the  matrix  elements  for  these  states. 

To  begin  with  we  shall  deal  with  the  BK  terms  for  the  triplet  and  singlet  1S2P^ 
and  1S2PX  states.  Although  the  procedure  of  applying  Fourier  transforms,  etc, 
is  identical  to  that  outlined  in  the  previous  section,  several,  different  types  of 
integrals  occur  which  are  worthy  of  mention.  These  integrals  occur  in  the  same 


( t-  *')  +  r  y  ( l  -*'1  +  ‘ili-y )  .1 


T  able 
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The  following  is  a.  listing  of  the  values  of  H 


CBk) 


for  capture  into  the  triplet 


1S2S  state  as  computed  from  eqn,  (5-1). 


g 

S  »  Zsf 

S\  3cxz 

0.  0 

-3. 2504 

-1, 3879 

=0, 6344 

=0, 3127 

=0, 0838 

0,  1 

-2, 8000 

-1. 0621 

=0,4243 

=0, 1780 

=0, 0277 

0,  2 

-1. 8076 

=0,4996 

=0, 1355 

=0, 0298 

0,0101 

0 ,  3 

-0, 8937 

-0. 1622 

=0, 0212 

0, 0060 

0, 0081 

0,4 

-0, 3436 

=0, 0447 

=0, 0009 

0, 0048 

0,  0031 

0 ,  5 

-0. 1022 

=0, 0166 

=0, 0008 

0, 0015 

0, 0010 

0,  6 

=0, 0241 

=0, 0103 

=0, 0014 

0, 0002 

0, 0002 

0,  7 

=0, 0073 

=0, 0072 

=0, 0013 

=0, 0002 

=0, 0001 

0,  8 

=0, 0057 

=0, 0046 

=0, 0009 

=0, 0003 

0,  9 

=0, 0050 

*0. 0027 

=0, 0007 

=0, 0003 

1,0 

=0, 0030 

=0, 0016 

=0, 0005 

=0, 0002 

1,  1 

=0, 0007 

=0, 0010 

=0, 0005 

1,2 

0, 0011 

=0, 0008 

=0, 0004 

1,3 

0,  0021 

=0, 0007 

1.4 

0, 0025 

=0, 0006 

1,5 

0, 0023 

=0, 0005 

1,6 

0, 0019 

1.7 

0, 0014 

1.8 

0, 0009 

1,9 

0, 0004 

2,  0 

0, 0001 
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T  able  4 


The  following  is  a.  listing  of  the  values  of  for  ca.pture  into  the 

triplet  1S2S  state  as  computed  from  eqn„  (5  =  2). 


g 

iS-  3c*4 

S\  4  c*2- 

0.  0 

=0. 1462 

=0. 1075 

=0. 0678 

=0. 0404 

0.  1 

=0. 1353 

=0. 0934 

=0. 0558 

=0. 0317 

0.  2 

=0. 1103 

=0. 0661 

=0. 0354 

=0. 0182 

0.  3 

=0. 0843 

=0. 0440 

=0. 0212 

0.4 

=0. 0643 

-0. 0299 

=0. 0129 

0.  5 

=0. 0504 

=0. 0207 

=0. 0079 

0.  6 

=0. 0406 

=0. 0144 

=0. 0047 

0.  7 

=0. 0330 

=0. 0098 

=0. 0027 

0.  8 

=0. 0267 

=0. 0065 

=0. 0014 

0.  9 

=0. 0213 

=0. 0043 

1.0 

=0. 0167 

=0. 0027 

1.1  -0.  0129 

1.2  “0. 0098 

1.3  =  0.  0074 

1.4  =0.0054 

1.5  =0.0039 
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The  ratio 


100 


is  plotted  against  the  energy  of 


the  incident  ion„ 


E  (mev.  ) 
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o 
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Figure  4 

The  cross  section  0”°=-  'sT^z.  +  Cn^^s  is  illustrated  by  the  solid 

curve  For  comparison  the  cross  section  0^.  (or  ) 

is  shown  by  the  dashed  curve  __  _  ___  __  °  The  experimental  points  of 

Allison  are  indicated,, 


-32- 


0.3 


0.6  — E(mev. ) 


1.5 


1.8 
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fashion  as  did  the  integral  lj/  of  eqn.  (4-'5)0  The  general  form  of:  the  integral 
which  occurs  is 


Ja  ~  -I  L^+Af+Y1]' 


(5  =  7) 


where  A  is  a  linear  combination  of  the  vectors  A  and  B  (i0  e0  =  ^c\  *■  ) 

and  ,  q  ,  and  are  defined  by  the  directions  of  A,  B  X  A,  and  B  respective¬ 
ly.  A  method  of  evaluating  these  integrals  is  outlined  in  Appendix  4. 

In  the  final  results  which  are  shown  below  the  upper  sign  (when  two  signs 
appear)  applies  to  capture  in  the  appropriate  triplet  state  and  the  lower  sign 
applies  to  the  singlet  state.  Thus  for  the  1S2P^  states  we  have 

H(  60  =  B  ['['['  yz 

7o  J o  Ja 


X  I  i 


be 


■*-  ± 


OlC. 


- ,  _  . 


Ck (  i-Ct  1  *i  cU  d-y  <1 
C/Sab 

whereas  for  the  singlet  and  triplet  1S2P^  — 


BK  matrix  elements  we  have 


A  -,4 


HUO  =  A  f  f  y2Ci~*)(A- % 

|\J  Jo  Jo  ^0 

bc» 


a  c. 


(5  =  8) 


w  i  ^  ^  '  _i_  -  —  b  (j — Lid — 

X  Up* s+'Uf1  t/W*^  +  C/WUV* 


„  a( 


^  ch  Y  <h 


(5-9) 


’•  -  \ 
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where  the  notation  in  both  (5-8)  and  (5-  9)  is  the  same  and  is  defined  as  follows; 


CX  =  XAS  )  b  -  z'  )  a.  Z'1  ^ 

where  the  appropriate  values  of  and  2-"  (i.  e,  whether  for  singlet  or  triplet 

states)  are  listed  in  Section  3, 


C  \  =•  *C\-*y)  ^  ; 

=  Azc,(\-C^  4-  , 

=  B*cvt\-0  4-  ^'C7_C\-C-l) 

iat  -  o?  %>iz  +  4-  b^U-y)  4- 

i  W  “  t^xyir  +  ^  (.i-x)  -4-  Q>  £  Li-y)  A  y 


K]  = 


,3/z.\  Ca-v-bV 


12  ^2ak)  L^t+M8j>y 


'/2 


(5=10) 

Again  the  evaluation  of  the  corresponding  nuclear  matrix  elements  is  straight” 
forward  but  tedious,  We  note  that  the  symbols  }  ^  ^  ^  and  ^3^  are 

still  defined  by  eqns,  (5  =  6),  (5  =  3),  (5  =  5),  and  (5=4)  respectively.  Also  to  differ- 
entiate  between  the  triplet  and  singlet  contributions  we  have  introduced  a  symbol 

rO 

botr  defined  such  that 

r—*  r 

6  -  ^  2  >  ^  4  for  singlet  terms 


and 


^  -  +1 


OJo  -  +  1 


=  -1 


L  -  1^2. 


t  *  3^  4 


for  triplet  terms, 


Thus  for  the  singlet  and  triplet.  1S2P-  states  we  have 


H^uc)  =  C  B 


/•  1 


6-1 


'o  J 


w( 


"f 


6  C : 


33C1 


dxdf  d; 


[  F*  F/fi*  /3,4  r*aP*  (5.ii 


Tl-T7.=  Ci  Cl-CtXl-x)(\-y)C\-?-)  5 

Th  --  Ta'-  > 


where 
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while  for  the  singlet  and  triplet  1S2P^  states  the  nuclear  matrix  elements  take  on 


-  A 


i 


£ 


r  ' 

"  \ 

Jq  - 

0  < 

0  F;3/*- 

!4Ci. 


where 


d x  c!^ 


T/'-Tj1  -  c.c/O-xXv-yrti-z-) 

t3'  -  v  -  yCi 


Now  only  the  numerical  integrations  on  a  computer  need  be  done  to  complete 
the  evaluation  of  the  cross  section  for  capture  into  the  singlet  and  triplet  1S2P 


states. 


. 
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-VI-  Possible  R  efinements 


The  results  of  the  last  two  chapters  were  obtained  by  an  application  of  the  first 
Born  approximation*  Obviously  the  use  of  the  second  Born  approximation,  which 
is  equivalent  to  writing  the  operator  U(  h)  ,  defined  by  eqn*  (2  =  5),  as 

U(M  |  -\-  Gr0('K\\J  (6=1) 


would  be  a  formidable  task.  Hence  one  is  led  to  search  for  procedures  which  are 
deemed  improvements  over  the  first  Born  approximation  but  which  are  still  amen  = 
able  to  feasible  calculations „ 

As  mentioned  in  the  Introduction  the  problem  of  single  electron  capture  by  fast 

protons  in  hydrogen  has  been  treated  by  taking  into  account  the  distortion  produced 

on  the  incoming  or  outgoing  waves  by  the  electron-nucieus  potential  of  the  bound 

system*  Below  we  illustrate  the  matrix  elements  which  would  have  to  be  evaluated 

were  this  type  of  refinement  to  be  used  in  treating  two  electron  capture  by  fast 

alpha  particles  in  helium,,  The  system  of  subscripts  used  is  still  that  illustrated 

21. 

in  Fig*  1*  By  employing  the  Chew=Goldber ger  identity  and  following  the  same 
procedure  of  McDowell^  we  may  write  eqn*  (2  =  8)  as 

fa+> 

—  \  1  -2,')  Qq.') 

+  terms  negligible  at  high  energies  j  (6-2) 


where  -  .  +• 

^  c,'j'  ~  I  4-  - 

E.  -  Wo-V/^  4-  L  & 

Remembering  that  \  at  high  energies,  i*  e» 

of  the  electron  wave  function  by  the  internuclear  potential, 


there  is  little  distortion 
we  have 


Htt  =<^tM\yA6  +  vMw„2Kwftr+uAf-i)0-  (6“3) 
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In  the  limit  of  very  high  energies  60^  '  and  eqn.  (6-3)  is  equivalent 

to  the  first  Born  approximation  as  we  have  used  it  throughout  this  thesis,.  It  is 
clear  that  the  computation  of  Ufc>a  from  eqn.  (6-3)  would  be  a,  tedious  task  to 
say  the  least.. 

Another  possible  refinement  might  consist  of  using  better  variational  approxi¬ 
mations  to  the  bound  states  of  the  helium  atom.  However  in  view  of  the  complexity 
of  the  results  which  we  have  obtained  from  two-parameter  wave  functions  it  is  not 
hard  to  see  that  the  addition  of  more  parameters  to  the  trial  wave  functions  would 
increase  the  complexity  of  the  calculations  enormously.  In  fact  one  soon  reaches 
the  conclusion  that  due  to  the  many-body  aspect  of  the  problem  under  consideration 
here  any  refinement  on  a  first  Born  approximation  calculation  might  increase  the 
complexity  of  already  arduous  computations  to  perhaps  an  unjustifiable  extent. 

At  any  rate  the  merits  of  any  one  refinement  will  be  unclear  until  the  experimental 
gap  is  filled  in. 
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-VII-  Conclusions  and  Discussion 

It  has  been  seen  that  the  effect  of  employing  an  improved  atomic  wave 
function  for  the  computation  of  the  cross  section  for  simultaneous  capture  of 

X 

two  electrons  by  fast  alpha  particles  in  helium  does  not  make  a  substantial 
difference  over  the  same  cross  section  computed  with  the  simpler  wave  function,, 
The  small  discrepancy  between  the  cross  section  of  GR  and  our  Os,^p  must 
be  interpreted  as  arising  from  the  various  errors  inherent  in  both  methods 
of  computation,  Thus,  for  the  above  type  of  calculation;,  it  seems  unreason¬ 
able  to  bother  using  a  more  refined  ground  state  wave  function. 

Unfortunately  the  evaluation  of  matrix  elements  for  capture  into  the  singlet 
and  triplet  1S2P^  and  1S2P^  states  would  have  been  too  time  consuming  with 

the  computational  facilities  which  were  available.  However  since  at  energies 

3 

near  0,6  Mev,  the  cross  sections  for  capture  into  the  triplet  1S2S  state  were 
found  to  be  about  12%  of  the  ground  state  value ,  it  would  not  be  unreasonable 
to  expect  a  contribution  of  the  order  of  50%  over  the  energy  range  from  0,25 
to  2  Mev,  to  the  total  capture  cross  section  from  capture  into  excited  states. 

It  should  be  recalled  that  in  the  energy  range  25  -  100  Kev,  JS  calculated  that 
capture  into  excited  states  accounts  for  one-third  of  the  total  cross  section 
for  electron  capture  by  protons  in  atomic  hydrogen.  It  is  expected  that  the 
numerical  calculations  for  capture  into  the  singlet  and  triplet  ISZP^  and 
lSZP^  states  will  be  completed  in  the  near  future. 

One  may  look  forward  eventually  when  reliable  experimental  cross  sections 
have  been  measured  and  the  excited  states  cross  sections  computations  are 
complete  to  a  settlement  of  the  controversy  over  when  the  nucleus-nucleus 
interaction  need  not  be  included  in  the  perturbation  Hamiltonian,  Certainly 
the  inclusion  of  this  term  as  in  the  protons  in  hydrogen  case  makes  a  large 
difference  in  the  calculated  cross  section. 


- 

. 

->  '  :  '■  t>a  /!MV 

.  c 


-39- 


Appendix  1 


A  general  method  of  performing  integrations  of  the  type  P„  Q9  R  and  S 

I  7  7 

listed  in  Section  II  is  to  replace  the  term  — - =rr  by  its  equivalent 

|T\  -  l\l 

1  1  ^ 


r  * — j  v  p 


\ 


^  ^  \X. 


z 


r,  z_j  v  r  )  Pjl  © ) 

1=0  1 


where  Q  =,  CUv$j  (,  f,  )  and  hence 

*  P^Cc^©,)  ^CC^02) 

^7  C-P-  1  pro  ^ 

-v  2^  Zj  (i  +■  l  Z  Cc^a,)  r  Cc^ coffin . 
\  X 

It  is  clear  that  if  in  integrals  of  the  form 

T  „  f  „R^»)  P^C  iA2.)  -V-  +_ 

1-  J  |r,  -rj  ^ 

either  F^  or  (  or  both  )  are  spherically  symmetric  functions  only  the 

jC  ~  0  terms  will  contribute  to  the  integral  and  hence  the  integrations  become 
trivial. 

The  integral  P  can  be  essentially  resolved  into  three  parts  denoted 

+'v\)  , 

e  1  tr  „ 


Ko  *■  J  IF.-F,. 

K  --  f—  ^ 

II'',-  hi 

It 2.  -  !  e^' cb,  3>z 


Using  the  above  method  K0  is  easily  evaluated  as 

(o<yZ  4- 


Kc  =  3ZTT2- 
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whence 


Kn"  (r  0 


n  B 


(  O 


Integrals  occuring  in  the  evaluation  of  Q  are  of  the  form 


Q 


m  n 


t  Q 


,Qri  p  0^)2. 


These  are  readily  evaluated  by  considering 


I 


I  _  /  c  r,"  +- 

Q^-J  17,1^1  e 


r  >  ^  2) 


whence  QJj(3  =  Kc  is  used  to  yield  the  result 

Q^rv  =  ("0 


a*r ^ 


Q 


oo 


The  validity  of  this  procedure  is  easily  checked  by  evaluating  a  few  of  the 
by  an  alternate  method  „  i0e.  by  expressing  the  terms  occuring  in 
the  integrand  by  their  respective  Fourier  transforms  and  then  integrating  the 
resultant  expression. 

The  integral  R  s  apart  from  a  factor  introduced  by  the  cm3  <9r  term  is 
identical  in  form  to  „  In  S  however  only  the  Jl-  I  term  in  the  expansion 


of 


I  ?,  -  M 


will  contribute  and  the  integral 


h,  -- 


In  -  ^ I 

is  readily  evaluated  yielding 


-ckI'T-vO 

£  CM  6)  Qf,  a>t 


- 


3  * 
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Appeadix  2 


The  following  three  dimensional  Fourier  transforms  were  found  useful 
in  reducing  the  BK  and  nuclear  matrix  elements; 
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Appendix  3 


All  forms  of  the  Feynman  auxiliary  integral  may  be  obtained  from  the 
simple  identity 

dx 


r 


ex  b 


[clx  -b  L(i-x)  J 


Repeated  differentiation  with  respect  to  a  and  b  gives  the  general  form 


1 
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N  (\-xV~  dx 


By  writing 


a  be 


=.  I  m  b  n  -  \  V _ 

tvn-OU^"')  !  JD[(Xx  -v  bOx)]^ 

dx  r'  dx 


and  then  using  the  integral  representation  for 

l  r  I 


o 

zc  we  get 
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-  zf  f  r 

JaJo  \r.  v.,  ,  1 


Again  repeated  differentiation  yields 
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Similarly  one  obtains  the  useful  representation 
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Appendix  4 


The  integral  defined  by  eqn„(5-7)  is  readily  evaluated  by  employing 


jJ 


the  integration  over  the  azimuthal  angle  will  vanish  and  we  are  left  with 

}  "  CM'c  BjAN> 

=  -tt’-Cj.  B  U-t-O- 

(  JD-  \)  \  ^  2.4. '3 

With  the  aid  of  the  coordinate  system  below  in  which  /\  lies  in  the  =  O 
plane  and  is  the  azimuthal  angle  of  ^  in  the  plane  perpendicular  to 

A  we  may  express  as 


%y 


, 

'■  '  '  •  ■  *  '  ....  :•  i 


-44- 


i  -rcY) 

Hence  due  to  the  integration  over  Q'  the  integral  J  '  is  zero,,  The 
projection  ^  may  be  written  as 

ir\j 


V  "  S 


—  9*  (  Aj  &  )  iw,  ~C  ^  +•  (  A)  &  )  T" 


Thus 
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